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An extended f/(3)L U{?>)r chiral effective field theory which includes pseudoscalar and vector 
meson nonets as dynamic variables is presented. The theory combines a hidden local symmetry ap- 
proach with a general procedure of including the rj' meson into chiral theory, and accounts for direct 
and indirect symmetry breaking effects via a mechanism based on the quark mass matrix. The 
theory is applied to anomalous radiative decays of pseudoscalar and vector mesons, using particle 
mixing schemes corresponding to different symmetry breaking assumptions and uniquely determined 
by the lagrangian presumed. Radiative decays of light flavor mesons are best explained within the 
framework of a one mixing angle scheme and provide evidence for SU{3)f and nonet symmetry 
breaking. 
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^ ' I. INTRODUCTION 

(N : 

' The decays of light flavor mesons have been discussed vigorously in the literature |2^. Particularly, phenomeno- 
^ ■ logical models based on effective field theories have been rather successful to explain anomalous and non-anomalous 
T^lj- ] processes involving vector and pseudoscalar mesons, tensor and higher-spin mesons, J/tp decays into a vector and a 
(N| ■ pseudoscalar meson, and many other related decays and topics [p| jHl^ , [l^ , |l5t . Particularly successful are effective field 
^-H ' theories based on the Hidden Local Symmetry (HLS) approach [p5| , claiming to provide a consistent framework which 
. explains a vast amount of data, and more importantly, reflects on the energy properties of the more fundamental 
QCD lagrangian. 

^ ', In the chiral limit, the QCD lagrangian exhibits an SU (3)l (S* SU{3)ii symmetry which breaks down spontaneously 

a to SU{3)v, giving rise to an octet of pseudoscalar light Goldstone bosons. The QCD spectrum contains a ninth 
singlet boson because the axial U{1) symmetry is broken by the anomaly. Nowadays it is well accepted that the 
I lightest pseudoscalar mesons (7r+,7r~, tt", K+, K~, K°, K° and ??) are candidates of the Goldstone octet. Though 
Qh! considerably heavier, the 77'(957 MeV) meson is considered to be the pseudoscalar singlet. 

^ • Our main interest in the present work is to develop and apply a J7(3)l C^(3)_r chiral theory where the the 
pseudoscalar and vector meson nonets play the role of dynamic degrees of freedom. The theory combines the HLS 
approach of Bando et al. with a general procedure of including the 77' meson into a chiral theory p6|~p9t, and is 
constructed in close analogy with QCD. The SU{3)l ^SU{3)ii local symmetry based QCD lagrangian is extended 
by adding a term proportional to the topological charge operator (the so called winding number density), which gives 
■ rise to the well known Wess-Zumino-Witten (WZW) anomaly term and, renders the lagrangian to be U{3)l ^ U{3)r 
locally symmetric p^ . Most importantly, the effective lagrangian constructed, exhibits the fundamental symmetries of 
this extended QCD lagrangian and accounts for direct and indirect symmetry breaking effects via a mechanism based 
on the quark mass matrix. A natural way to account for SU{5) flavor symmetry breaking terms involves an expansion 
of the QCD generating functional in powers of the quark masses. Here as well we use a similar expansion to generate 
symmetry breaking terms. We believe that the theory proposed defines an accurate framework for electroweak and 
strong interactions of light flavor mesons and can be extended easily to include tensor and higher spin mesons as 
well as baryons. Although, the structure of the lagrangian seems similar to the original model of Bando, Kugo, 
Yamawaki |2^] and variants of this model |^,|l^,^, it includes additional terms due to including the 77' and the 
mesons as dynamical variables. Namely, the trace of the nonet field matrices are nonvanishing and therefore, the 
symmetric as well as the symmetry breaking lagrangian parts must include the square of the traces of vector and 
axial-vector covariants, in addition to terms proportional to traces of covariants square. These additional terms play 
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an important role in determining the particle mixing schemes and the meson mass spectra. In marked difference 
with previous studies [^| Jl^ , [l^ the meson masses, mixing angles and other observables like radiative decay widths or 
form factors are calculated straightforwardly from the lagrangian presumed without any additional assumptions. The 
theory guarantees reproducing the meson masses and as a result the mixing angles are functions of the symmetry 
breaking scales and not free model parameters. 

Our paper is organized as follows. In section II we construct the lagrangian. In order to be fully consistent with 
QCD, a symmetry breaking companion is added to each of the unbroken lagrangian terms. These symmetry breaking 
companions affect the kinetic and mass terms of the lagrangian which, in turn determine the mixing schemes of 
state and coupling constants |Q. Two alternative ways are proposed, corresponding to SU{S)f symmetry breaking 
(Alternative I) and U(3) p symmetry breaking [Alternative II). The strong condition of nonet symmetry often 
used to relate the octet and singlet pseudoscalar states is not presumed a priori. In the limit of nonet 

symmetry (in the sense of equal singlet and octet radiative decay constants, Fq = Fg) the Alternative II scheme 
reduces to a scheme equivalent to that of the quark flavor basis (QFB) of Feldmann et al. jl0-12|. The realization 
of these mixing schemes is described in section III. In section IV we apply our model to study anomalous processes. 
Namely, we calculate radiative decay widths of y° — > P'^j, P° V^j and P° — > 77, with — 7T,K,ri,r]' and 
= p, K*,uj, 4>. Global fit to data is performed in order to determine numerically the symmetry breaking scales and 
evaluate the success of the different schemes to explain data. We summarize and conclude in section V. 



II. THE EFFECTIVE LAGRANGIAN 



The main objective in the present section is to construct a most general C/(3)l ® C^(3)_r local symmetry based 
effective lagrangian, for pseudoscalar and vector meson nonets interacting with external electroweak fields. As a 
non-linear representation of the pseudoscalar nonet fields we define 12§,ET 




[/(P,ryo + Po^?)=r(P,r;o + Fo^9) = exp<^i^P + i\/-— ryol , (2.1) 



where riQ{x) stands for the pseudoscalar singlet and P{x) for the pseudoscalar Goldstone octet matrix, 

(2.2) 
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with obvious notation. In this representation, the octet degrees of freedom are contained in the unimodular part of 
the field U while the phase detU = expi(V6??o/-Fb) involves the singlet only. The vacuum angle serves here as an 
auxiliary field ^(x) which renders the variable X{x) = —ilndetU + 1} = ^/6t]o{x) / Fq + •d{x) to be invariant under 
U{3)l C^(3)_R transformations p6| , p7t . Since there exists no dimension-nine representation for U{3)l ^ U{3)r, the 
octet (Fg) and singlet (Fq) radiative decay constants may well be different. We follow Refs. pl[-p3[| and define vector 
and axial-vector covariants, 

^M = ^[^^5^e] + ^(^^+^^^^0 ' (2-3) 

A. = ^{Ct,9,C} + ^(eV^e-e^.e^) • (2.4) 

Here and l^ represent the standard model external gauge fields; = -I- and = Vfi — o,fn with and 
being the vector and axial vector external electroweak fields, respectively. For pure electromagnetic interactions 
Ifj. = Tfj. — —eQAfj_ where denotes the electromagnetic field and Q — diag{2/3,—l/3,—l/3) is the quark charge 
operator. 



Under U (3) l <S> ^ (i) r the field, Eqn. 2.1, transforms as, 

U' = P[/Lt, (2.5) 

with R e [/(3)fl, L £ U{3)l. The vector (F^) and axial- vector (A^) covariants transform, respectively, as gauge 
and matter fields, i.e., 

F;, = KT^K^ + iKd^K^ , (2.6) 

a;, = K^^K^ , (2.7) 
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where K(U, R, L) is a compensatory field representing an element of the conserved vector subgroup ?7(3)y |3l|-|33|]. 
The dynamical gauge bosons are defined as a 3 x 3 vector field matrix which transforms as, 

= KV^K^ + '-Kd^K^ . (2.8) 

Clearly, the vector — gV^ as well as the axial vector transform homogeneously, and to lowest order (i.e. 
with the smallest number of derivatives) the lagrangian can be constructed from the traces of , (F^ — gV^)^, 
^Ai: i^i-i — gVfj.), D'^'d, and arbitrary functions of the variable X{x) all of which being invariant under U{3)l ^ U{3)ii 
transformations. Then to lowest order, a most general form of a symmetric effective chiral lagrangian is. 
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with, 



L = LA + aLv- -TriV^.V") , (2.9) 



La = WiiX)Tr{A^A'^) + Wi{X)TriA^)TriAn + 

W5{X)Tr{A^)D^'d + We{X)D^-dD^d , (2.10) 
Lv - Wi{X)Tr{[T^ - 5^^]^^ - .g^) + 

Wi{X)Tr{T^ - gV^)Tr{T^^ - gV^') , (2.11) 

and, 

D^d = d^-d + Tr{r^-l^) , (2.12) 
V^, ^ d^V, - d,V^ - ig[V^, K] . (2.13) 



All three terms of the lagrangian L in Eqns.2.9-2.11 are invariant under U(3)l ^ U{3)r transformations. 

The coefficient functions, Wi, and Wi, are constrained by parity conservation to be even functions of the variable 
X. In addition by requiring that the normalization of the HLS kinetic term be equal 1/2, it is easy to show that 
Wi{0) = Fg, Wii^O) = {F§ — Fg)/S, and We{0) = l/2. Such a normalization ensures that the pseudoscalar singlet 
couples to the singlet axial current with a strength Fq while the octet states couple to the octet axial currents with a 
strength Fg. 

Although the lagrangian, Eqn. |2.9| , appears similar to that of Bando et al. j2^, the terms La and Ly are differ- 
ent. Namely, including the 77' and (j) mesons as dynamical variables involves additional terms with Tr(A^)Tr(A^), 
Tr{A^)D^§ and Tr(r'' - gV^')TT{T^ - gV^). Clearly, such terms must be included in any J7(3)l ® C/(3)fl symmetry 
based effective field theory, where the traces of A^ and — gV^ do not vanish. In addition we have introduced the 
coefficient functions Wi{X) which becomes constants in the SU{i) symmetry limit. The kinetic term of the pseu- 
doscalar mesons as well as their strong and electroweak interactions with the Goldstone fields are all included in La- 
As for the vector mesons, Ly incorporates all interactions with the pseudoscalar fields. The kinetic term is written 
explicitly as —^Tr{Vp^yV'^'^). Similar to the lagrangian of Bando et al. |^^, the sum La + aLy contains, amongst 
other contributions, a vector meson mass term ^ V^V^ , a vector-photon conversion factor ~ V^A^^ and the coupling 
of pseudoscalar pairs to both vectors and photons. The latter coupling can be eliminated by choosing a value a — 2, 
which allows incorporating the conventional vector-dominance in the electromagnetic form-factors of pseudoscalar 
mesons [^, and eliminating the coupling of a pseudoscalar particle to two photons. Data analysis ||l^ with variants 
of BKY lagrangian argue for a value a = 2.4. This parameter is determined anew |34| using our lagrangian. 



A. SU{3)f symmetry breaking 



As already indicated in the introduction, the expansion of the QCD generating functional in powers of quark mass 
term results with SU(3) flavor symmetry breaking terms. There is no unique way to introduce such terms into the 
lagrangian, but it is quite natural choosing them to be (i) Hermitian, (ii) proportional to powers of the quark mass 
matrix and, (iii) recover the unbroken lagrangian smoothly in the limit of vanishing symmetry breaking parameters. 
To be consistent with QCD, we add a symmetry breaking companions to each of the unbroken lagrangian terms, 
using the a procedure as described below. First, the Goldstone meson mass degeneracy is removed by adding a 
U{3)l (S) U{i)ii symmetry violating mass term. Most generally such a term reads, p^ , |27| , p9[| , 

Lrn = -Wo{X) + W2{X)Trx+ + iW3{X)Trx- , (2.14) 
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with, 



X± = 2Boimk ± ^^M^^) , Bo = ml/im^ + m^) , (2.15) 

where mu,md,ms are the up down and strange quark masses, and A4 ~ diag{mu,md,ms). Next, in order to 
incorporate symmetry breaking corresponding to each of the unbroken La and Ly lagrangian densities, we define a 
universal Hermitian matrix B, 

B = -x+ ■ (2.16) 

4:Bo{mu + md + ms) 

Then, symmetry breaking companions for La are constructed in two alternative ways. The first (hereafter referred 
to as Alternative L) breaks octet symmetry {SU{i)F) only, and corresponds to a quadratic form of the Goldstone 
meson kinetic energy term. Let, 

C/g = = exp(*^P) (2.17) 

be the pure octet field matrix and let, 

\^\[ild^^s} + \(ilr^is-isUl) , (2.18) 

be the octet axial vector covariant. Then a general symmetry breaking lagrangian La would be. 

La = Wi{X) {cATr{B A^A^) + dATr{BA^BA^)) + 

Wi{X)dATr{BA^)Tr{BAi') + W5{X)cATr{BA^)D>'§ . (2.19) 

The second alternative {Alternative II) breaks U{2>)f symmetry and uses the nonet axial vector A'' ( instead of 
as above). Namely, 

La = Wi{X) {cATr{B A^Af^) + d^rr(SA^BA'^)) 
+W4X) {cATr{BA^)Tr{A^') + dATr{BA^)Tr{BA>')) + 

W5(,X)cATr{BA^)D^'§ . (2.20) 

This expression gives a bilinear meson kinetic energy term. Similarly the asymmetric companion of Ly would be, 

Lv = W^{X) {evTr{B[T^ - gV^][T^ ~ gV'']) + 
dvTriB[T^ - gV^]B[T'' - gV''])) + 
WiiX) (cvTriTf, - gV^)Tr{B[r^ ~ gV^]) 

+ dvTriB[T^~gV^,])TriB[T^-gV^])) . (2.21) 

In the expressions above, CA,dA, cy and dy are symmetry breaking scales to be determined from data analyses. It is 
to be stressed that La and Ly differ also from the ones defined in Refs. ||25|,|5|,|l^ . First, like our symmetric La and 
Lv, the asymmetric companions La and Ly involve additional terms which are absent in the SU(3) limit. Secondly, 
in the original work of Bando et al. and variants of this model jlj], the scales di are taken to be equal to (ci/2)^, 
based on their assumption that the physical fields are related to the bare fields through rescaling, i.e., 

^ph = . (2.22) 

Contrary to this assumption, Bramon et al. have assumed di = 0. In the our model, the field rescaling is dictated 
by the diagonalization procedure and is uniquely determined by the form of the lagrangian ]30|] . For the more general 
lagrangian constructed above, and in order to reproduce the actual spectra of the physical mesons the constants 
Ci and di must then be treated as independent parameters. Thirdly, unlike Bando et al. [p5| the symmetry breaking 
matrix B is Hermitian. It is similar (but not identical) to that of Bramon et al. j|] and rather close to the one proposed 
by Benayoun and O'Connell flj]. Albeit, the lagrangian is constructed in close analogy with QCD and allows for 



symmetry breaking in a universal manner. Particularly, the matrix B, Eqn. 2.16 , enables us to maintain the QCD 
ratios of isospin to SU{3) symmetry breaking scales. 
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Summing all terms, including the well known Wess-Zumino-Witten term Lwzw |35 3^, the lagrangian assumes 
the form, 

L = La + La + L,n + a{Lv + Ly) - ^Tr{V^,V'''') + Lwzw + ■■■ , (2.23) 



where " . . ." stands for terms accounting for the regularization of one loop contributions [I3|,[4y2q] . Note that Ly includes 
vector meson mass terms which also violate symmetry. Generally speaking, a symmetry breaking companion for the 
vector meson kinetic term —jTr{Vi_i^V'^'^) should have been added also. However, at present, no evidence exists for 
such asymmetric terms and therefore will not be considered in the following discussion. The lagrangian constructed 
above should describe the mass splitting of the pseudoscalar and vector mesons (see for example p7[ ). Since the 
ground state of the field U is proportional to the unit matrix, one can set the auxiliary field 1} — pM. With these 
simplifying assumptions the quantities Wi and Wi become functions of the singlet field (770) only. To lowest order 
their expansions read. 

Wo = const + Fi (^wo^ + .. ^ , (2.24) 
W^=Fi(^^ + w^^ + ..)j , (2.25) 
= f + + , (2.26) 

W^3 = ^(^3^ + ...), (2.27) 
W, = ^l^(^l + w,^ + ..}j (2.28) 

(2.29) 
(2.30) 
(2.31) 
(2.32) 

Here the coefficients Wi, Wi and Wi are free parameters not yet determined. In the applications to be discussed below, 
however, only a few combinations (rather than each) of these parameters are needed. 

III. SYMMETRY BREAKING AND PSEUDOSCALAR MESON MIXING 

Symmetry breaking terms affects the form of the kinetic and mass lagrangian densities of the Goldstone mesons 
which in turn uniquely determine the state and decay const ants mixing schemes pO]]. In what foll ows w e realize the 
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mixing schemes, corresponding to the Alternative I ( Eqn. 2.19 ) and the Alternative II ( Eqn. 2.20) lagrangians. 
These two alternatives allow us to evaluate SU{3)f versus U{3)f symmetry breaking. A third scheme corresponding 
to the QFB of Feldmann et al. [0, is also realized from the Alternative II by requiring nonet symmetry {Fq = Fg). 
Generally, for any EFT the kinetic and mass terms can be written in the form, 

ifcm = i(a^$)/C(9^$) + ^<fMH , (3.1) 

where $ represents the intrinsic meson field matrix, /C and are the kinetic and mass matrices. In the presence of 
symmetry breaking IC and are non-diagonal, and Eqn. 3.1 does not have the standard quadratic form as invoked 



by the Klein-Gordon equation for the physical fields. It can be reduced into this standard form by applying three 
consecutive steps, which transform the intrinsic fields into the physical meson fields according to ]30| ], 

$ = e%h ; e = TRn . (3.2) 



5 



Here T represents a unitary transformation which diagonaUze the kinetic matrix /C, R stands for rescaling of the 
fields, and Q is another unitary matrix which diagonahze the resulting mass matrix. Let us work in details a general 
scheme for the mixing of two states, where presumably, the i] and 77' are linear comb inations of an octet (rjs) and a 
singlet (rjo), only. Using the symmetry breaking Alternative II lagrangian, Eqn.|2.20|, the kinetic and mass terms are 
bilinear and the corresponding matrices read. 



Here, 
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(3.3) 
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mj(2m + TTig) 
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(3.4) 



and r — F^/ Fq is a measure of nonet symmetry breaking. For simplicity , small terms of order ~ CATfi/ms, dAtn/ms 
are neglected and isospin symmetry (m„ — rrid = m) is assumed in the expressions listed above. Note that /C and 
have non-diagonal 2x2 submatrices, and as demonstrated in Ref. |3^ the transformation Q involves two mixing 
angles and two rescaling parameters. First as in Ref. [3C| ] we diagonalize K, using the unitary transformation. 



K 
Vs 
Vo 



== T 




T = 



10 
10 

cos A sin A 

— sin A cos A 



(3.5) 



This leads to. 



with, 



'^0 = 7: 



H88 + noo - \J (koo - i^s&Y + 4k|| 

K88 + KOO + \l (kqo - Kss)^ + 4k|| 



and, 



tan A 



'^80 



1 — tan" A Koo — Kga 



(3.6) 

(3.7) 
(3.8) 

(3.9) 
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Next we rescale the fields (tt, K, rjg,, ?7o) into (tt, K, r)^, t)q) using, 



R = diag{\, zk, z, f) = diag{\, 



1 



1 



1 



(3.10) 



/kk ^/l^8 ^/Ko 

Following these two steps, the kinetic term acquires the standard quadratic form and the resulting mass matrix is. 



RT-^M^TR 



ill 
ill 



ills 

/iio M§o 



with, 



~2 2 

~2 2 2 



— z"^ (cos^ -^l^ss + ^i^^ -^Moo ~ 2 sin A cos A/^gp) , 
Aso = ^/ (sinAcos A(^|8 - fi^o) + (cos^ A - sin^ A)^|o) , 
Aoo ~ (sin^ '^Mss + -^Moo + ^ ^™ '^'^^ ^/^so) • 
As a last step we diagonalize the mass matrix, Eqn. ^.11 , by applying a second unitary transformation, 



cos X sm X 
— sinx cosx 



and take the eigenvalues to be equal to the masses of the physical particles, i.e., 



m-l^ fil , 

2 _ 2 2 



and. 



tanx 



^88 Moo 



1 — tan X 

By definition the physical 77 and 77' fields are orthogonal and are related to the intrinsic fields via. 



where. 



e' 



z COS A cos X — / sin A sin x z cos A sin x + / sin A cos x 
— z sin A cos x ~ / cos A sin x —2: sin A sin x + / cos A cos x 



(3.11) 



(3.12) 
(3.13) 
(3.14) 
(3.15) 
(3.16) 



(3.17) 

(3.18) 
(3.19) 

(3.20) 
(3.21) 

(3.22) 

(3.23) 

(3.24) 



Using this transformation, we may rewrite the nonlinear representation, Eqn. 2.1 in terms of the physical fields, e.g., 

/2 

C/ = expi^7', (3.25) 
where V stands for the pseudoscalar nonet matrix, 
^ + -ig(X„ry + X,,V) 



TT 



7r+ ZsK+ 



(3.26) 
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The coefficients Xi and li, (i = rj^rj') are listed in Table |[ It is worthy to mention that by taking the eigenvalues of 
the resulting mass matrix to be equal to the masses of the physical mesons we have eliminated explicit dependence 
on the model free parameters Wi. Consequently, the mixing angles A and x, the rescaling parameters z and /, the 
coefficients Xi and Yi^ and the pseudoscalar nonet matrix V become functions of just three symmetry breaking scales 
CA, cLa and r. 

This enables us to incorporate symmetry breaking indirectly in the pseudoscalar nonet matrix, Eqn. 3.26| . We refer 



Where, 



to these as indirect as opposed to direct symmetry breaking due to the broken La and Ly lagrangian companions. 

At this stage some comments are in order : 
(i) A better known expression of the transformation Q is , 

= cos^ X + sin^ X; z cos x = zi cos , / sinx = — zi sin V'l, (3.28) 

= cos^ X + -2^^ sin^ Xi / cos X = Z2 cos 7/12 , zsinx = — Z2sin7/;2, (3.29) 

/ P 
tan -01 = tanx, tan -f/;! = ^ tan '02 , (3.30) 

z 

e,, = A - 02 , - A - 01 . (3.31) 

This exact general expression of the field transformation involves two mixing angles On and 9^' as well as two field 
rescaling parameters zi and 2:2, and is dictated by the kinetic and mass terms of the lagrangian presumed. A similar 
two angle scheme was proposed by Escribano and Frere ||l8|,^ , however their transformation does not account for the 
field rescaling and violates orthogonality of the r] and 77' states. 



z. 



(ii) With the Alternative I symmetry breaking lagrangian La, 2.19, the kinetic term assumes a quadratic form, 
with Kgg = 1 + 2{ca + 'iy!i)/3, kqs = 0, and kqq = 1. AH other matrix elements of K. and M"^ remain the same as in 



Eqn. 3.4. Consequently, z = l/-\/l + 2{ca + dA)/'3> and / = 1, and the angle A vanishes, so that the Alternative I 
scheme involves a single mixing angle. 

(iii) In the limit of exact nonet symmetry, r — 1 [Fs = Fq), the Alternative II pseudoscalar meson kinetic matrix 
becomes diagonal in the so called quark flavor basis (QFB), where the flavorless mesons are represented by fields with 
or without strange quark content (77s ^ ss and r]q ~ {uu + dd) / . With this basis the kinetic matrix is diagonal, 
which leads to a one mixing angle scheme. Note that Feldmann et al. |p^-|l^ do not account for field rescaling and 
therefore, their mixing angle is equivalent to our mixing angle X: which defines the transformation of the rescaled 
fields into the physical ones. 

(iv) As a result from the diagonalization procedure, in all mixing schemes, the rescaling parameters and mixing 
angles are functions of the symmetry breaking scales and can not be taken to be free parameters, as commonly 
assumed in previous studies p^jr^ , p^ . 

To conclude this section, we briefly consider the vector meson mixing problem. Here, the kinetic energy term has 
the standard quadratic form. Thus one needs diagonalizing the mass matrix only, which effectively depends on four 
parameters nriy = 2F^ag, cy, dy and w^, where my stands for the vector meson nonet symmetric mass. All four 
parameters can be fixed by taking the calculated physical meson masses to be equal to their experimental values. 
Following a similar diagonalization procedure as above, it is straightforward to show that the physical vector nonet 
matrix with non-ideal mixing reads, 

^{p^ + uj + e^) p+ K*+ 

p- ^[-p^ + u; + e^) K*^ I . (3.32) 

K*- K*° 0-61 

with e = 0.047 ± 0.001, being a measure of the non-strange (strange) admixture in (w). Nearly the same value of e 
is deduced from the width of the — > 7r7 decay (see next section) . In practice, the vector meson mixing angle turns 
to be very close to the one of ideal mixing, and likewise, the corresponding strange (non-strange) admixt ures i n the 
nonet vector matrix are small. Since we worked in the exact isospin symmetry limit, the vector field, Eqn. 3.32| , does 



not account for p — lo mixing. The derivation of the broken matrix is straightforward and will be given elsewhere jMl 
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IV. RADIATIVE DECAY WIDTHS 



The main objective in the foUowing section is to apply our theory to calculate decay width of the anomalous 
processes listed in tabl e fS/\ Details of the calculations are described in previouspublications and will not be repeated 
( see for example Ref. |15|). Our analyses generalize the treatments of Refs. |^,|,^jl5|, in the sense that we treat all 
nonet mesons on equal footing and most importantly we use mixing schemes which are uniquely determined by the 
lagrangians presumed. 

In our analyses we use pseudoscalar and vector field matrices expressed in terms of physical fields and therefore 
allow for indirect and direct symmetry breaking effects. The lagrangian is factorized in the form, 

Lpyy = Lpy^ + C]^Lpj^, (4.1) 
Lyp-r — Lypj + CwLyp^ , (4.2) 

where Lp^-y {Lyp~f) represents the overall contribution of the unbroken anomalous lagrangian to the P77 (VP^) 
interaction, and Lp^-y (Lvp-y) is the corresponding direct symmetry breaking companion. Here cw stands for the 
direct symmetry breaking parameter. As mentioned above, both Lp^^ and Lpyy account for symmetry breaking via 
the pseudoscalar and vector meson matrices. To write these terms explicitly, we start from the anomalous lagrangian 

Lanomalous — Lyyp + Lwzw{Pll)- (4.3) 

In terms of the covariants A^, — .9V'^, Viiu and = dfjTn — di,T^ — i[r^,ri,] (see sect II for definitions), the 
pseudoscalar- vector- vector (PVV) coupling, Lyyp, at lowest order has at most six contributions. 



L'ylp = gie'"'"''Tr(V^^,[F„ - ^r^]Ap) + g^e^^-^'^^TriT^^iV^ - ^r„]A^) + 
gse''-'-^Tr{V^,.)Tr{[Vo. _ ir„] A^) + 54e^''"^rr(r^,)Tr([F„ - ir„]A/3) + 
g,e''-'-^Tr{V^,.)TriVa, „ ir„)Tr(A^) + 

^g,P-/3r^(r^^)Tr(K. - ir,)Tr(A^) , (4.4) 

where gi, {i = 1,...6) are arbitrary coefficient functions of the variable X. By rearranging terms we may write 
Lanomalous in a compact form, e.g., 

LvP^ = gv^e'"'''l^d^A,Tr{Q{do,Vp, V}) , (4.5) 

Lp^^ ^ gp-^e^'''''^d^A,daAf,Tri{Q^,V}) . (4.6) 

Here gp and gv represent certain combinations of the function coefficients gi and constants in Lwzw- As in section 

II, the direct symmetry breaking terms are constructed by inserting the quantity B, Eqn. p. 16 , in the expressions 
above, i.e., 

Lvp^ - gv^e'''"'^d^A,Tr{Q{B,{do^Vp,r}}) , (4.7) 

-Ttv 

Lp-ry = gp^^e^'-'^d^A.daApTriiQ'^iB^r}}) . (4.8) 



A. The V — > P7 and P Vj Processes 

The relevant vertices for a vector (pseudoscalar) meson decaying into a pseudoscalar (vector) meson and a photon 
are, 

ViVP^) = -igv^viVP)e'^-''^''k,ei'^^p^ej,''^ , (4.9) 
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where elf^ (p) and el^^ (k) are the polarization (four-momentum) of the vector meson and final photon, respectively. 
The quantities v incorporate all internal symmetries of the processes under discussion and are listed in Table [|. In 
terms of these vertices the widths of the decays V Pj and P V'f are, 

r{vp^)^Gv^-^^^^^^HVPr , (4.10) 



riPV^)^3Gy^-^^^^^^^\v{VPr , (4.11) 



with, 



B. The P 77 decays 



The vertices for decays of a pseudoscalar meson into two photons are. 



e 



2 



V(P77) = -2igp^viP)e'^''"''k,^e^jh2J^^ , (4.13) 



Fs 



with el7'' and Ca"^ being the polarizations of the final photons, and fci, k2 their corresponding four-momenta. Again, 
the functions v{P) contain the internal symmetry information and are listed in Table |l[ With these vertices the 
decay rate is given by. 



r(P77) = Gp^|z;(P)P, (4.14) 



where. 



TT / e 



2\' ^nr.^2 



C. Numerical Analysis and Results 

The expressions quoted above for the decay widths involve three symmetry breaking scales {ca, Aa, r), a direct 
symmetry breaking scale (cw), and two coupling constants (f;y, gp). We fix the coupling constants from experiment 
and we treat the other quantities as free model parameters to be determined from global fit. As a general rule the 
experimental masses and decay widths used are the best fit values reported in the latest review of particle properties 
IpTj . The pion radiative constants is taken to be = Fg = 93MeV. 

Consider now the values of the coupling constants. From the experimental width of the uj decay, and Eqn. 



4.10| , one obtains r(u;7r7) = Gv{ml - mlf/{mlFi) = (716 ± 43)^6^ This relation yields, 

Gy = (1.44 ±0.04) . 10"^ , gy = 0.220 ±0.006 . (4.16) 

In practice, these same constants are obtained from a similar relation, r(p7r7) = G{m^—m'^)^/ (m^F^) = (76±10)iireT^ 
p7| for the p ^ ttj decay. Similarly, from the decay width of tt ^ 77 one finds, r(7r°77) = 9Gpm^/2F| = 
(7.8±0.55)ey 

Gp = (4.91 ±0.07) • 10"^ gp = 0.073 ±0.001 . (4.17) 

Furthermore, the vector meson strange-nonstrange admixture parameter e can be deduced from the decay widths of 
the and m mes ons into a pion and a photon. The^— > ttj decay is forbidden in the limit of ideal mixing. However, 



from Eqns. 4.10 and the experimental decay rates |37[| one obtains. 
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{ml - ml)m^ 



H 3 



(5.8 ± 0.6)is:ey 
(716 ± 43)if eT/ 



= 0.008 ±0.001 



(4.18) 



This ratio gives e = 0.043 ±0.004 , a value which is equal within experimental error, to that deduced from vector mass 
matrix diagonalization(see section II). The calculated decay widths and parameter values as obtained from global fit 
to data using different mixing schemes are summarized in Tables III and [V. Based on the / dof values listed, the 
one mixing angle Alternative I scheme provides the best explanation of the data. This observation remains valid, 
should we have used an older set of data |^^, though the resulting fit qualities are slightly poorer. 

Alternative I -One mixing angle scheme First we recall that in this case f — 1, kqs = and A = 0. This 



simplifies Eqns.3.5 3.10,3.17, 3.24 and allows expressing the mixing angle in terms of the physical meson masses. 



tan dp 



l-tan^gp ^ 



2M; 



- 1 



■88" 



where /ilgz^ is determined by the pion mass, kaon mass and rescaling parameters, i.e. 



- 1 



(4.19) 



(4.20) 



The rescaling parameter z must not be considerably smaller than zk- Indeed, from Eqn. 3.13 for the kaon mass and 
Eqns. 3.21 for the 77 and r/' masses we may write. 



4m|. 1 

3 7Th-f^ ^ 



1 ^ 1 
3 m,^ ~ 



(4.21) 



or equivalen tly, 1.03z > zk. Thus z is at most a factor of ~ 3% smaller than zk- With the Alternative I parameter 
set of Table III, the rescaling parameters and mixing angle are, 



ZK = 0.87 ±0.04, 



0.89 ±0.06, 6ip = -(14.3 ±2.2)° 



(4.22) 



How significant are the departure of these parameters from their values at the exact limit of U{3)l^U{3)u 
symmetry? Clearly, the exact SU{S)f symmetry limit values cw = 0, ca = cIa = 0, r = 1 would be inconsistent 
with the measured value of the ratio r(K*'^ K'^j) /r{K*~^ K~^j); with cw = this ratio becomes 4 as opposed to 
the experimental value of 2.34 ± 0.43. In addition, direct symmetry breaking alone is not sufficient; with ca = 
and Cw = —0.2, the calculated width T{K*'^ K'^j) is about 30% higher than experimental value, well beyond the 
measurement accuracy. Based on fit quality both direct (i.e. cw 7^ 0) and indirect (i.e. CA,dA 0,r ^ 1) symmetry 
breaking terms are needed to explain data. 

As a further check we may estimate these scales using data directly. First, from the ratio, 



= 4 



1 



jCw 



I - Cw 



(117± 10)ii:ey 
{50±5)KeV 



= 2.34 ±0.43 



(4.23) 



one obtains cw = — 0.19±0.04 in full agreement with the global fit results. Next, the decay width r{K*''^K'^'-f) involves 
the kaon rescaling parameter (or equivalently the scales ca)- One finds, zk = 0.86 ± 0.08 and ca = 0.52 ± 0.22. The 
other two parameters are deduced from the ratios, T{(j)ri'y)/T{ujr]j) and, T{(j)ri^) /T{ri' p^), which yield, 



= -0.45 ±0.2 , r = 0.98 ±0.1 



(4.24) 



With this set of parameters the corresponding mixing angle Op ~ —(16.2 ± 2.4)° and /dof — 36/6 is more than a 
factor of 10 higher as compared to the Alternative I fit results. A global fit to data with the assumption of an exact 
nonet symmetry, i.e. r = 1, gives cw = -0.18 ± 0.05, ca = 0.64 ± 0.06, = -0.32 ± 0.04 and Op ^ -(15.2 ± 2.2)° 
but with x^/dof = 13.2/7. Again the fit quality is significantly poorer as compared to the Alternative I results. 

It is worthy of mention that the mixing angle is very sensitive to the ratio of the rescaling parameters z and 
Zk- The smallest mixing angle is obtained with z = zk (i.e., c?a = — ca/4) and readily grows for increasing z/zk- 
However, a global fit with these parameters taken to be equal is again far poorer with x^/dof = 8.2/7 and cw = 
-0.26 ±0.05, CA = 0.56 ±0.06, r = 0.86 ±0.1, 9p = -(8.2 ±2.0)°. The mixing angle is far less sensitive to the direct 
symmetry breaking scale cw- Based on these global fit analyses we may conclude that nonet symmetry breaking and 
field rescaling parameters depart slightly but significantly from unity. 
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Finally, by comparing Eqn. 4.20, with the familiar quark model result 

,2 / „^ 2 



the symmetry breaking measure is, 



"^i8 = ^(4^-l))(l + A), (4.25) 



A - 1 . (4.26) 



The Alternative I parameters correspond to A = 0.05 ± 0.01. 

The Alternative II -Two mixing angle scheme The mixing angles and rescaling parameters corresponding 



to the Alternative II fit parameters of Table [II arc 



A = (25.1 ± 2)°, X = -(34.9 ± 3)°, zk = 0.95 ± 0.06, z = 1.00 ± 0.06, / = 0.88 ± 0.06 , or equivalently, 

61,, = -(5.8 ±2)°, 61,,. = -(12.7 ±2)°, zi = 0.98 ± 0.05, Z2 = 0.96 ± 0.05 . (4.27) 

In view of the very poor fit quality (confidence level of less than 0.05), it seems quite unjustified to use a two mixing 
angle scheme in analyzing radiative decays. 

Quark flavor basis scheme As indicated already, this scheme corresponds to the Alternative II in the limit 
of exact nonet symmetry {r — \). In this case the values of ca and (Ia are considerably different, and the rescaling 
parameters and mixing angle are, 

(40.0 ± 2.0)°, = 0.76 ±0.04, = 0.88 ± 0.05 . (4.28) 

With the standard octet-singlet mixing angle defined as 9p — (p ~ Oideai as in Ref. ||l^ we have 9p = —(14.7 ± 2.0)°. 
This value agrees remarkably well with the Alternative I results. Yet, with jdoj — 19.8/7 (confidence level less 
than 0.01) this scheme like the Alternative II is far poorer than the Alternative I. 

It is of interest to note that similar mixing angle values were extracted by several authors using various phenomeno- 
logical models. Here we mention few examples. Feldmann et al. ||l^,|l^ using their QFB analysis reported a value 
Op = —(14.8 ± 2.9)°. Cao and Signal obtained value of 9p = —(14.5 ± 2.0)° from analyzing large momentum 
e+e" ri,ri' ^ 27. Somewhat a less negative value, 6p = — (11.59 ± 0.76)°, was determined by Benayoun et al. 
|E5|,E6| . The fact that different model analyses predict similar results should not be surprising since the mixing angle 



is sensitive to ratios of scaling parameters (see Eqn. 4.19) rather than to their actual values. In fact, since this ratio 



is close to one, the mixing angle should be determined rather accurately from the experimental meson masses. This 
also means that in similar analyses where rescaling is neglected or the corresponding ratios of parameters are close to 
unity, as in our Alternative I case, similar values of 9p are expected. However as our analyses proves, rescaling plays 
an important role in explaining the decay widths. 



V. SUMMARY AND DISCUSSION 



In this paper, we have constructed a local C/(3)i U{3)p symmetry based effective field theory, using a generalized 
version of the HLS approach of Bando et al. ||2^ and a general procedure of including the rj' into a chiral theory [p6|-p9| . 
At lowest order the lagrangian comprises two terms La-, Ly describing respectively, the interactions of pseudoscalar 
and vector meson nonets, and a "kinetic" term for the vector mesons. In variance with the BKY model |2^ and 
variants of this model these terms and likewise their asymmetric companions La and Ly involve terms which 

are proportional to the square of traces, in addition to terms proportional to the trace of vector square and axial- vector 
square. Such terms must be included in any U{3)l^ U{3)p symmetry based theory, and play an important role in 
reproducing the actual meson masses. In a way similar to that proposed by Bramon et al. M and Benayoun et al. 
1^,^, we have constructed "direct" asymmetric companions La and Ly using a matrix B which is proportional to 
quark mass matrix. With this choice of B, our theory predicts the same ratios of isospin to SU (3) symmetry breaking 
scales as in QCD. Our lagrangian includes also the exact form of symmetry violating mass term as derived by Herera- 
Siklody et al. | |29[ |. Two alternative forms of symmetry breaking parts were proposed, the first of Alternative I breaks 
SU{3)f only while the second Alternative II breaks U{3)f- Each of these alternatives corresponds to a unique state 
mixing scheme, which can be derived through the diagonalization procedure of the corresponding kinetic and mass 
terms, a procedure which turns to be very useful in constraining the theory parameters. A third scheme corresponding 
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to the QFB of Feldmann et al. , is realized from the Alternative II by requiring nonet symmetry. As a result of 
the diagonalization of the kinetic and mass terms of the lagrangian, the mixing schemes, i.e. the mixing angles and 
field rescaling parameters, depend on three symmetry breaking scales (c^,fiA,?') only, and by all means can not be 
treated as free parameters. This relation between the mixing schemes and the lagrangian form has been overlooked 
in previous studies |l^Jl^ , p^ . 

Using these three schemes we have applied our theory to the pseudoscalar sector. The case of the vector meson 
sector is considered elsewhere . The results obtained reflect upon the nature of the symmetry breaking required 
to explain radiative decay widths. Based on global fit analyses, we observe that the Alternative I provides by far 
a better description of the data considered. This argues for S'[/(3)f symmetry breaking supplemented with broken 
nonet symmetry. The other alternatives of ?7(3)f symmetry breaking with or without nonet symmetry breaking yield 
rather poor fits and seem unjustified. 

There have been numerous data analyses in the last three decades attempting to deduce a reliable and accurate value 
of the pseudoscalar mixing angle [p]J^J^ P, p^ , p7| , pO[ . The values reported range from 9 — —23° to as high as 6 = —10°. 
In marked difference with these previous studies, in the present work the particle mixing scheme is well related to the 
lagrangians presumed. All observables are calculated from the lagrangian with no additional assumptions. Albeit, we 
believe that the theory proposed provides an accurate framework for the study of electroweak and strong interactions 
amongst light flavor mesons. 
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Alternative I 
z cos 9p ~ \/2r sin Op 

-2z cos 9p — \f2r sin 6'p 

2 sin 6'p + \f2r cos Bp 

-2z sin + \/2r cos 6p 



Alternative II 



QFB 

\/3 cos (j) 



X'r, 



Y' 



z cos A cos x~ f sin A sin x+ 
\/2r(— zsinAcosx — /cos Asinx) , 

— 2{z cos A cos X ~ .f sin A sin x) + 
\/2r{—z sin A cos x ^ / cos A sin x) , 

z cos A sin x + / sin A cos x+ 
\/2r(— 2 sin A sin X + / cos Acos x) , 

— 2(z cos A sin x + / sin A cos x) + 
\/2r(— 2 sin A sin x + / cos A cos x) ■ 



— ^&/ks si 
\/3 sin ( 



•y/ Q/ks cos (j) 



TABLE I. The Xi and Yi coefficients, Eqn. 3.26. Mixing angles are 6p for the Alternative I and 
Ks = 1 + ca + dA is a matrix element of the kinetic matrix in the QFB. 



for the QFB scheme. 



v{pTl) = i 

vliuit) — 1 



\zk{1 + ^cw) 



v{K*+K+) = v{K*-K-) = i2K(l - cw) 



Xri 



v{u>r]) 



)(r,) = ^[5X, + (l+2cw')K,] 



v{Lori ) 



y{H)--^Yr,>(l + cw) 



viv') = 7S [5X,y + (1 + 2cw)Y^' 



TABLE IL The internal symmetry factors v(VP) and v(P), Eqns. 4.9, 4.13 





Cw 


CA 


dA 


r 


xVdof 


Alternative I 


-(0.20 ±0.05) 


(0.64 ±0.06) 


-0.25 ±0.04 


0.91 ±0.04 


3.1/6 


Alternative II 


-(0.27 ±0.05) 


0.2 ±0.05 


0.1 ±0.02 


0.94 ±0.05 


18.6/6 


QFB 


-(0.19 ±0.05) 


(1.4 ±0.1) 


-1.1 ±0.1 


*1 


19.8/7 



TABLE IIL Symmetry breaking scales and x^ jdoj from global fit to data. Values marked with an asterisk were kept fixed. 
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Decay 


rea:p(KcV) 


Alternative I 


rcaic(KeV) 
Alternative II 


QFB 


O — ^ TTT 


76 ± 10 


76 ± 10 


76 ± 10 


76 ± 10 


LO — > TT'y 


716 ± 43 


* 716 


* 716 


* 716 


r 11 


36 ± 12 


40.0 ± 4.1 


47.4 ±4.2 


52.8 ± 4.5 


uj — Ti'y 


5.5 ± 0.85 


5.2 ± 0.45 


6.1 ± 0.5 


6.8 ± 0.6 


Y '1 I 


57.8 ± 1.6 


60.8 ± 2.4 


65.7 ± 3.5 


60 ± 3 




+0 20 
n 30 ^"-^^ 

-0.16 


0.37 ±0.03 


0.24 ±0.03 


0.4 ±0.05 


T] ~* p-y 


61.2 ± 7.5 


70.7 ±5.9 


70.3 ±4.1 


78.5 ±6.6 


rj' Lwy 


6.12 ±0.75 


6.47 ±0.45 


6.4 ±0.4 


7.2 ±0.5 


7r" 77 


(7.8 ±0.55)10"^ 


* 7.8 • 10"^ 


* 7.8 ■ 10"^ 


* 7.8 • 10"^ 


?7 ^ 77 


0.460 ± 0.050 


0.468 ± 0.03 


0.58 ±0.02 


0.637 ±0.03 


rj' ^ 77 


4.280 ± 0.280 


4.02 ± 0.3 


3.68 ±0.24 


4.52 ±0.3 




117± 10 


106 ±8.3 


120.7 ±4.6 


113 ±8 




50 ±5 


48.6 ±3.2 


63.3 ±5.2 


48 ±5 


xVdof 




3.1/6 


18.8/6 


19.8/7 



TABLE IV. Calculated radiative decay widths. The Alternative I, Alternative II and QFB correspond to the parameter 
sets of Table The widths marked with an asterisk were used to evaluate the coupling constants gp and gv (see text). The 
experimental decay widths of the second column are the best fit values of Ref. . 
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